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SUMMARY 

The  characteristic  properties  and  evolutionary  processes  of  naturally  occurring  and  man-made  aerosols  is  a  sub¬ 
ject  of  considerable  importance  for  modeling  electro-optical  propagation  through  the  atmosphere.  The  Navy’s  effort 
is  centered  upon  predicting  propagation  conditions  that  impact  on  sensor  and  communication  systems  in  the  marine 
environment.  Since  the  aerosol  properties  of  an  air  mass  depend  primarily  upon  its  source  region  and  upon  the  alter¬ 
ation  it  undergoes  in  transit  to  the  observation  site,  effective  methods  for  air  mass  tracing  are  essential  for  modeling 
the  atmospheric-optical  environment.  A  series  of  five  electro-optical  field  programs  conducted  at  San  Nicolas  Island 
(1978-1980)  afforded  an  opportunity  to  compare  two  independent  methods  for  identifying  air  mass  regimes.  One 
technique  utilized  measurements  of  atmospheric  radon  concentrations  as  an  air  mass  tracer  while  the  other  approach 
employed  a  computerized  objective  analysis  scheme  to  automatically  construct  surface  air  trajectories.  The  categori¬ 
zation  of  over  270  trajectories  into  low  (<  10  pCi/m3),  medium  (>  10,  <30  pCi/m3),  and  high  (>30  pCi/m3)  radon 
concentrations  provides  an  approximate  description  of  the  air  parcel  history  for  these  conditions. 

The  lowest  radon  concentrations  occurred  when  the  air  traversed  a  northwesterly  to  westerly  oceanic  path 
before  ending  at  San  Nicolas  Island.  The  highest  radon  levels  were  associated  with  air  parcels  that  either  passed  over 
land  or  through  the  southern  California  bight.  In  a  few  cases  relatively  high  radon  concentrations  were  observed  for 
over-water  trajectories,  indicating  earlier  overland  travel  or  horizontal/vertical  mixing  with  continental  air.  The 
synoptic  features  and  meteorological  factors  responsible  for  the  observed  results  are  discussed. 
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1.  INTRODUCTION 


San  Nicolas  Island  (SNI)  is  located  off  the  southern  California  coast  about  120  Km  west-southwest  of  Los 
Angeles  (figure  1)  and  is  routinely  used  by  the  Navy’s  Pacific  Missile  Test  Center  as  a  weapons  testing  and  radar 
tracking  site.  In  a  concerted  effort  to  characterize  the  effects  of  a  maritime  atmosphere  on  optical  and  infrared  trans¬ 
mission  pertinent  to  weapons  RDT&E,  the  U.S.  Navy  has  established  an  atmospheric  transmissometer  and  micro- 
meteorological  facility  on  the  northwestern  end  of  SNI  (figure  2).  Measurements  began  in  1978  under  the  sponsor¬ 
ship  of  the  Navy’s  Electro-Optical  Meteorology  (EOMET)  and  Optical  Signatures  (OSP)  programs.  To  date  a  total  of 
five  Atmospheric-Optical  measurement  programs  have  been  conducted  on  SNI  (Blanc,  1978* ,  19792,  1 98 1 3 , 19824; 
Matthews  et  al.,  19795;Jeck,  19796;  Rosenthal  et  al.,  19797; Clark et  al„  1981a8, 1981b9, 1982a10, 1982b11). 

One  of  the  major  objectives  of  the  EOMET  program  is  to  characterize,  for  a  variety  of  atmospheric  conditions, 
the  marine  boundary  layer  meteorological  parameters  that  influence  electro-optical  propagation.  One  approach  in 
this  endeavor  is  to  derive  meteorologically  dependent  aerosol  models  (e.g.,  Bamhardt  and  Streete,  197012;  Wells  et 
al.,  197713;  Katz  et  al.,  197914)  from  which  extinction  coefficients  can  be  computed  using  Mie  scattering  theory. 
Since  the  refractive  index  of  aerosol  particles,  which  is  required  in  Mie  computations  of  aerosol  extinction,  depends 
primarily  on  air  mass  history  (Richter  and  Hughes,  1981 15),  it  is  reasonable  to  attribute  at  least  some  of  the  marine 
aerosol  model  deficiencies  encountered  in  validation  efforts  (Hughes  and  Richter,  197916;  Noonkester,  198017; 


1  Blanc,  T.  V.,  1978:  Micrometeorological  data  for  the  Cooperative  Experiment  for  West  Coast  Oceanography  and  Meteorology 
(CEWCOM-78)  at  San  Nicolas  Island,  California  Vol.  I  and  II.  NRL  Memorandum  Report  3871,  Aug.  1978.  UNCLASSIFIED 

2Blanc,  T.  V.,  1979:  Micrometeorological  data  report  for  the  November  1978  Electro-Optics  Meteorology  (EOMET)  experiment  at 
San  Nicolas  Island,  California.  NRL  Memorandum  Report  4056,  Aug.  1979.  UNCLASSIFIED 

3Blanc,  T.  V.,  1981:  Report  and  analysis  of  the  May  1979  marine  surface  layer  micrometeorological  experiment  at  San  Nicolas 
Island,  California.  NRL  Report  8363,  Dec.  1981.  UNCLASSIFIED 

4 Blanc,  T.  V.,  1982:  The  data  base  for  the  May  1979  marine  surface  layer  micrometeorological  experiment  at  San  Nicolas  Island, 
California.  NRL  Report4713,  May  1982.  UNCLASSIFIED 

^Matthews,  G.  B.,  B.  E.  Williams,  A.  Akkerman,  J.  Rosenthal,  and  R.  de  Violini,  1978:  Atmospheric  transmission  and  supporting 
meteorology  in  the  marine  environment  at  San  Nicolas  Island -semiannual  report.  PACMISTESTCEN  Technical  Report  TP-79-19, 
Dec.  1978.  UNCLASSIFIED 

^Jeck,  R.  K„  1979:  Aerosol  particle  size  measurements  at  San  Nicolas  Island  during  CEWCOM-78.  NRL  Memorandum  Report 
3931,  Mar.  1979.  UNCLASSIFIED 

7Rosenthal,  J.,  T.  E.  Battalino,  V.  R.  Noonkester,  and  H.  Hendon,  1979:  Marine/continental  history  of  aerosols  at  San  Nicolas 
Island  during  CEWCOM-78  and  OSP  III.  PACMISTESTCEN  Technical  Publication  TP-79-33,  Apr.  1979.  UNCLASSIFIED 

8Clark,  R„  T.  E.  Battalino,  R.  Helvey,  and  J.  Rosenthal,  1981a:  Atmospheric  conditions  and  airmass  history  at  San  Nicolas 
Island  during  the  November  1978  EOMET/OSP  experiment.  PACMISTESTCEN  Technical  Publication  TP-81-30,  Dec.  1981. 
UNCLASSIFIED 

9Clark,  R„  T.  E.  Battalino,  and  R.  Helvey,  1981b:  Atmospheric  conditions  and  airmass  history  at  San  Nicolas  Island  during  the  30 
April-11  May  1979  EOMET/OSP  experiment.  PACMISTESTCEN  Technical  PublicationTP-81-33,  Dec.  1981.  UNCLASSIFIED 

l0Clark,  R„  T.  E.  Battalino,  and  R.  Helvey,  1982a:  Atmospheric  conditions  and  airmass  history  at  San  Nicolas  Island  during  the 
14-25  January  1980  EOMET/OSP  experiment.  PACMISTESTCEN  Technical  Publication  TP-82-21,  May  1982.  UNCLASSIFIED 

1 'Clark,  R.,  T.  E.  Battalino,  and  R.  Helvey,  1982b:  Atmospheric  conditions  and  airmass  history  at  San  Nicolas  Island  during  the 
22-31  July  1980  EOMET/OSP  experiment.  PACMISTESTCEN  Technical  Publication  TP-82-22,  May  1982.  UNCLASSIFIED 

12Bamhart,  E.  A.,  and  J.  L.  Streete,  1970:  A  method  for  predicting  atmospheric  aerosol  scattering  coefficients  in  the  infrared. 
Applied  Optics,  9,  1337-1344.  UNCLASSIFIED 

13Weils,  W.  C.,  G.  Gal,  and  M.  W.  Munn,  1977:  Aerosol  distributions  in  maritime  air  and  predicted  scattering  coefficients  in  the  infra¬ 
red.  Appl.  Opt.,  16,654-659.  UNCLASSIFIED 

14Katz,  B.  S.,  K.  Hepfer,  and  N.  E.  MacMeekin,  1979:  Electro-optics  meteorological  sensitivity  study.  NSWC  TR-79-67,  Apr.  1979. 
UNCLASSIFIED 

'^Richter,  J.  H.,  and  H.  G.  Hughes,  1981:  Electro-optical  atmospheric  transmission  effort  in  the  marine  environment.  NOSC  Tech¬ 
nical  Report  696,  May  1981.  UNCLASSIFIED 

16 Hughes,  H.  G.,  and  J.  H.  Richter,  1979:  Extinction  coefficients  calculated  from  aerosol  size  distributions  measured  in  a  marine 
environment.  Proceedings  of  the  Society  of  Photo-Optical  Instrumentation  Engineers,  San  Diego,  CA.,  Aug.  29-30,  1979. 
UNCLASSIFIED 

1 7Noonkester,  V.  R„  1980:  Offshore  aerosol  spectra  and  humidity  relations  near  southern  California.  Second  Conference  on  Coastal 
Meteorology,  Los  Angeles,  CA.,  Jan.  30-Feb.  1,  1980.  UNCLASSIFIED 
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Figure  2.  Sen  Nicolas  Island  and  the  Marine  Environment  Test  Range  (METR). 


Hughes,  198018)  to  air  mass  differences.  The  next  generation  of  marine  aerosol  models,  such  as  the  “Navy  Model” 
(Gathman,  198219),  will  explicitly  incorporate  air  mass  effects. 

Two  independent  methods  for  characterizing  air  mass  regimes  were  utilized  during  the  EOMET/OSP  measure¬ 
ments  at  SN1.  One  technique  constructed  surface  air  trajectories  for  parcels  terminating  at  SN1,  and  the  other 
approach  measured  the  concentration  of  atmospheric  radon  gas  for  use  as  an  air  mass  tracer.  The  purpose  of  this 
report  is  to  compare  these  two  air  mass  identification  techniques  and  to  discuss  some  of  the  related  meteorological 
factors. 


2.  ATMOSPHERIC  RADON 


Radon  222,  a  product  of  the  uranium  238  decay  series,  is  an  inert  radioactive  gas  (half  life  3.82  days)  pro¬ 
duced  by  the  decay  of  radium  226  (half  life  1622  years)  which  is  found  in  various  quantities  in  the  so  md  rocks  of 
the  earth’s  crust.  Most  of  the  rare  gas  emanates  from  continental  areas,  with  less  than  2%originat>.  i  die  oceans 
(Wilkening  and  Clements,  197520).  After  formation  in  the  ground,  radon  diffuses  into  the  atmo  <-e  where  its 
short  lived  daughter  products  become  attached  to  aerosol  particles  that  are  transported  with  the  p  Jmg  winds. 
The  majority  of  the  radon  burden  collected  by  an  air  mass  depends  upon  its  duration  over  land  ns  and  the 
radon  emissions  rate  from  those  regions;  the  removal  processes  are  governed  by  radioactive  decay,  di  and  tur¬ 

bulent  mixing  (Larson  et  al.,  1979c21).  Because  radon  is  produced  mainly  on  land  masses  and  has  t  •»-  .ong  half 
life,  it  has  been  widely  used  as  an  air  mass  tracer  over  oceanic  areas  (Fontan  et  al.,  196322;  Vilenskii  et  al.,  196723; 
Rama,  197024;  Prospero  and  Carlson,  197025;  Wilkniss  et  al.,  197426;  Moore  et  al.,  197427;  Subramanian  et  al., 
197728;  Larson,  1978  29,  1979b30;  Larson  et  al.,  197931;  Larson  and  Bressan,  198032;  Larson  and  Jeck,  198133). 


18Hughes,  H.  G.,  1980:  Aerosol  extinction  coefficient  variations  with  altitude  at  3.75  4m  in  a  coastal  marine  environment.  J.  Appl. 
Meteor.,  19,  803-808.  UNCLASSIFIED 

19Gathman,  S.  G.,  1982:  Navy  aerosol  model.  Paper  presented  at  the  Annual  Review  Conference  on  Atmospheric  Transmission 
Models,  AFGL,  Hanscom  AFB,  Mass.,  May  18-19,  1982.  UNCLASSIFIED 

20Wilkening,  M.  H.  and  W.  E.  Clements,  1975:  Radon  222  from  the  ocean  surface.  J.  Geophs.  Res.,  80,  3828-3830.  UNCLASSIFIED 

21Larson,  R.  E.,  D.  J.  Bressan,  K.  W.  Marlow,  T.  A.  Wojciechowski,  and  J.  L.  Heffter,  1979c:  A  comparison  of  concentrations  of  fis¬ 
sion  products.  Radon  222,  and  cloud  condensation  nuclei  over  the  North  Atlantic.  Pageoph.,  117,  874-882.  UNCLASSIFIED 

22Fontan,  J.,  D.  Blanc,  and  A.  Bouville,  1963 :  Meteorological  conditions  dependance  of  radon  concentrations  in  air  above  the  Atlan¬ 
tic.  Nature,  197,  583-584.  UNCLASSIFIED 

23Vilenskii,  V.  D.,  G.  V.  Dmitrieva,  and  Yu.  Kirasnopovtrev,  1967:  Natural  and  artificial  radioactivity  of  the  atmosphere  over  the 
ocean  and  its  relation  to  meteorological  factors.  U.S.  Atomic  Energy  Commission  Report  AEC-TR-6711.  UNCLASSIFIED 
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25Prospero,  J.  M.,  and  T.  N.  Carlson,  1970:  Radon-222  in  the  North  Atlantic  Trade  Winds:  its  relationship  to  dust  transport  from 
Africa.,  Science,  167,  974-977.  UNCLASSIFIED 

20Wilkniss,  P.  E.,  R.  E.  Larson,  D.  J.  Bressan,  and  J.  Steranka,  1974:  Atmospheric  and  continental  dust  near  the  antarctic  and  their  cor¬ 
relation  with  air  mass  trajectories  computed  from  Nimbus  5  satellite  photographs.  J.  Appl.  Meteor.,  13,  512-515.  UNCLASSIFIED 

27Moore,  H.  E.,  S.  E.  Poet,  E.  A.  Martell,  and  M.  H.  Wilkening,  1974:  Origin  of  222  Rn  and  its  long-lived  daughters  in  air  over  Hawaii. 
J.  Geophys.  Res.,  79,  5019-5024.  UNCLASSIFIED 

28Subramanian,  S.  K.,  C.  Rangarajan,  S.  Gopalakrishnan,  and  C.  D.  Eapen,  1977:  Radon  daughter’s  radioactivity  levels  over  the  Ara¬ 
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29Larson,  R.  E.,  1978:  Radon  222  measurements  during  marine  fog  events  off  Nova  Scotia.  J.  Geophys.  Res.,  83,  415-418. 
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Relative  changes  in  radon  concentration  are  more  meaningful  as  an  atmospheric  tracer  than  the  absolute  concentra¬ 
tion  (Larson  et  al.,  197931).  This  is  true  because  the  factors  affecting  radon  exhalation  depend  upon  the  variable 
conditions  of  soil  moisture  (Pearson  and  Jones,  1 96534),  atmospheric  pressure  (Clements  and  Wilkening,  I97435). 
crustal  radium  deposits  (Larson  and  Jeck,  198133),  and  the  strength  of  the  radium  source  (Wilkening.  197436). 
Radon  levels  of  a  few  pCi/m3  (picocuries  per  cubic  meter)  are  typical  over  oceanic  air,  while  concentrations  of  sev¬ 
eral  tens  of  pCi/m3  are  found  over  land.  An  increase  in  radon  concentration  over  ocean  areas  indicates  the  presence 
of  an  air  mass  with  a  recent  overland  history.  Thus  monitoring  fluctuations  in  radon  activity  provides  a  technique  for 
determining  whether  a  given  location  is  experiencing  an  intrusion  of  continental  or  maritime  air. 

The  radon  measurements  were  made  by  the  Naval  Research  Laboratory  at  the  micrometeorological  tower 
facility  which  is  located  on  the  northwestern  end  of  SNI  about  20m  from  the  high  tide  waterline  (Blanc,  1 978 1 ). 
The  instrumentation  employed  was  developed  by  Larson  (197337)  and  utilizes  a  measurement  technique  that 
assumes  equilibrium  between  radon  and  its  daughter  products  (polonium  218,  lead  214,  bismuth  214).  Radon  levels 
are  determined  by  collecting  the  daughter  products  on  filter  paper  and  measuring  the  beta  activity  with  plastic  scin¬ 
tillators  optically  coupled  to  photomultiplier  tubes  (Larson,  197337;  Larson  and  Bressan,  197838).  Although  some 
questions  have  been  raised  concerning  the  assumption  of  equilibrium  in  radon  daughter  products  (Shapiro  et  al., 
197839),  it  is  well  established  that  disequilibrium  in  maritime  environments  is  unlikely  except  during  rain  when 
daughter  products  can  be  selectively  removed  (Larson,  1979a40;  Gat,  196641). 


3.  TRAJECTORY  ANALYSIS  TECHNIQUE 

A  computerized  objective  mesoscale  analysis  technique  developed  by  Barnes  (197342)  has  been  adapted  by 
one  of  us  (R.  A.  Helvey)  to  provide  automatic  determination  of  air  parcel  trajectories  based  on  hourly  surface  wind 
observations  at  stations  distributed  over  the  region  of  interest.  A  trajectory  is  constructed  from  a  series  of  vectors 
representing  successive  hourly  displacements  of  a  particular  air  parcel.  Each  of  the  individual  displacements  is 
obtained  from  wind  components  computed  for  the  current  position  of  the  air  parcel  by  use  of  the  modified  Barnes 
technique. 

This  technique  interpolates  wind  observations  from  a  number  of  stations  to  specified  analysis  Joints  utilizing 
a  Gaussian  weight  function  of  the  form 


Wj  =  exp  (-4Rj2/Ds2) 

where  Wj  is  the  weight  applied  to  the  i  th  observation  for  the  analysis  point,  Rj  is  the  distance  between  the  i  th 
observation  point  and  the  analysis  point,  and  Ds  is  a  parameter  that  determines  the  response  characteristics  of  the 
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5290.  UNCLASSIEIED 
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Science,  183,413-415.  UNCLASSIEIED 
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weight  function.  Employing  the  exponential  space  weighting  function  allows  all  wind  observations  to  objectively 
influence  the  values  assigned  at  the  analysis  points.  The  analysis  algorithm  proceeds  in  two  steps.  First,  a  rough  esti¬ 
mate  of  the  analyzed  wind  is  obtained  using  a  relatively  large  value  for  Ds  which  effectively  minimizes  the  contribu¬ 
tions  from  small-scale  irregularities  in  the  observed  data.  Estimates  are  likewise  obtained  at  the  observation  points 
themselves;  these  estimates  are  then  subtracted  from  the  original  data  to  provide  a  correction  field.  In  the  next  step 
the  correction  field  is  interpolated  using  a  decreased  Ds,  which  now  permits  smaller-scale  features  to  influence  the 
analysis.  By  reducing  Ds  in  this  manner,  Barnes  found  that  satisfactory  pattern  resolution  could  be  achieved  with 
considerably  less  computational  effort  than  required  by  the  widely-used  Cressman  (195943)  technique.  For  the  sta¬ 
tion  distribution  in  our  study,  an  initial  value  of  60  nm  for  Ds,  reduced  to  15  nm  for  the  second  pass,  appeared  to 
yield  the  most  satisfactory  results.  The  final  analysis  value,  P,  for  n  observations  is  computed  from 

p  =  2  /  2  wj 

i=l  i=l 

where  Gj  is  the  weight  the  i  th  correction  receives  at  the  analysis  point  (Brandes,  197644). 

The  analyses  are  carried  out  on  a  Hewlett-Packard  9845T  Desktop  Computer,  which  has  a  CRT  (cathode-ray 
tube)  and  other  features  that  make  it  particularly  suited  for  interactive  data  entry,  error  detection  and  correction, 
and  display  of  intermediate  results.  A  number  of  steps  necessarily  precede  the  preparation  of  trajectories,  including 
entering  and  storing  station  identifiers  and  locations  and  hourly  wind  data  for  the  periods  of  interest.  Hourly  wind 
analyses  can  be  generated  to  assist  in  checking  for  errors  and  in  visualizing  wind  flow  patterns;  these  are  displayed  on 
the  CRT  in  the  form  of  wind  vectors  plotted  at  grid  points  over  the  domain  of  analysis.  Trajectories  are  plotted  at 
hourly  intervals  within  specified  periods;  they  can  be  determined  for  durations  up  to  72  hours  either  forward  or 
backward  in  time  for  specified  points  of  origin  or  termination.  Due  to  the  scarce  wind  information  available  in  our 
study  and  because  of  the  meteorological  uncertainties  (mixing  effects)  inherent  in  following  an  air  parcel,  the  trajec¬ 
tories  discussed  in  this  report  are  for  24-hour  duration.  Since  the  objective  was  to  determine  possible  effects  of  pre¬ 
vious  history  of  air  parcel  characteristics,  they  were  obtained  backward  in  time  from  a  given  termination  point  (SNI). 

Analyses  were  performed  for  a  330  x  260  nm  region  centered  near  SNI.  Hourly  surface  wind  directions  and 
speeds  for  12  southern  California  coastal  and  island  weather  stations  were  used.  In  addition,  winds  were  estimated 
from  surface  weather  charts  for  three  points  (PI,  P2,  P3)  over  the  open  ocean.  The  trajectory  computations  are  ter¬ 
minated  if  the  path  passes  outside  the  map  boundaries.  Topographical  and  air  mass  mixing  effects  are  ignored,  thus 
only  horizontally  advected,  near-surface,  air  motions  are  considered.  In  the  future,  low-level  cloud  motion  vectors 
derived  from  successive  geostationary  satellite  imagery  will  also  be  incorporated  to  obtain  estimates  of  winds  over 
data-sparse  oceanic  areas.  The  locations  of  all  1 5  sites  are  listed  in  Table  1 . 


4.  COMPARISON  OF  AIR  TRAJECTORIES  AND  RADON  LEVELS 

Figure  3a  illustrates,  for  the  five  EOMET/OSP  measurement  periods,  a  total  of  273  24-hour  trajectories  termi¬ 
nating  at  SNI  during  periods  when  radon  measurements  were  available.  Most  of  the  trajectories  originate  from  loca¬ 
tions  north  of  SNI  and  a  conspicuous  trajectory  void  exists  to  the  southeast.  Radon  amounts  were  observed  to  vary 
by  more  than  two  orders  of  magnitude;  the  largest  radon  level  recorded  was  126.1  pCi/m3  and  the  smallest  was  0.8 
pCi/m3.  The  large  variability  observed  in  air  parcel  history  and  radon  activity  (figure  3a),  which  results  primarily 
from  the  superimposition  of  different  air  mass  regimes,  can  be  resolved  into  trajectory  paths  that  correspond  to  dis¬ 
tinct  radon  concentration  categories.  Samson  (198045)  has  successfully  used  a  similar  approach  to  study  air  trajec¬ 
tories  and  sulfate  aerosols.  In  our  investigation  we  have  found  three  radon  concentration  categories  (low;  <10 


43Cressman,  G.  P.,  1959:  An  operational  objective  analysis  system.  Mon.  Wea.  Rev.,  87,  367-374.  UNCLASSIFIED 

44 Brandes.  K.  A.,  1976:  Hourly  surface  weather  maps  and  analyzed  fields  of  meteorological  variables.  NOAA  Technical  Memo¬ 
randum  ERL  NSSL-80,  Dec.  1976.  UNCLASSIFIED 

43Samson,  P.  J„  1980:  Trajectory  analysis  of  summertime  sulfate  concentrations  in  the  Northeastern  United  States.  J.  Appl.  Meteor., 
19.  1382-1394.  UNCLASSIFIED 
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Table  1.  Station  Location*  Utad  in  Trajectory  Computationi. 


Identifier 

Location 

Longitude 

(Watt) 

Latitude 

(North) 

SMX 

Santa  Maria 

120.417° 

34.933° 

SB  A 

Santa  Barbara 

119.833° 

34.433° 

LAX 

Lot  Angeles 

118.383° 

33.933° 

SNA 

Santa  Ana 

117.867° 

33.667° 

SAN 

San  Diego 

117.167° 

32.733° 

VBG 

Vandenberg  Air  Force  Bate 

120.567° 

34.733° 

NTD 

Point  Mugu 

119.117° 

34  177° 

NUC 

San  Clemente  Island 

118.583° 

33.017° 

NFG 

Oceanside 

117.367° 

33.300° 

SM4 

San  Miguel  Island 

120.363° 

34.033° 

SZ1 

Santa  Cruz  Island 

119.634° 

33.995° 

SN5 

San  Nicolas  Island  (Site  A) 

119.575° 

33.277° 

P1 

Open  Ocean 

122.000° 

35.000° 

P2 

Open  Ocean 

121.000° 

33.000° 

P3 

Open  Ocean 

119.000° 

32.000° 

pCi/m3;  medium:  >10,  <30  pCi/m3;  high:  >30  pCi/m3)  useful  in  distinguishing  between  the  air  flow  patterns 
depicted  in  figure  3a.  Figures  3b,  3c  and  3d  show  the  trajectory  paths  corresponding  to  the  low,  medium,  and  high 
radon  concentration  categories  respectively.  The  lowest  radon  levels  (figure  3b)  occurred  when  the  air  traversed  a 
northwesterly  to  westerly  oceanic  path  before  ending  at  SNI.  In  general,  the  highest  radon  concentrations  (figure  3d) 
were  associated  with  air  parcels  that  either  passed  over  land  or  through  the  southern  California  bight.  The  north¬ 
westerly  trajectories  are  observed  to  involve  brisk  regular  wind  speeds  while  other  directions  primarily  involve 
upstream  stagnation.  Trajectories  corresponding  to  the  medium  radon  category  (figure  3c)  consist  of  paths  that  are 
similar  to  those  found  in  both  the  high  and  low  radon  concentration  groups.  This  is  suggestive  of  a  transition  mode  if 
we  think  of  successive  changes  in  radon  reflecting  similar  changes  in  the  air  trajectories. 

Table  2  summarizes  the  trajectory/radon  sample  size  according  to  radon  category  for  the  individual  measure¬ 
ment  periods  and  figures  4,  5,  and  6  illustrate  the  trajectory  paths  for  each  period  and  radon  category.  In  figure  5a, 
relatively  high  radon  concentrations  were  observed  for  both  the  northwest  (parallel  to  the  coast  line)  and  southwest 
(towards  the  open  sea)  trajectories,  suggesting  earlier  overland  travel  or  horizontal/vertical  mixing  with  continental 
air.  The  two  overwater  trajectories  passing  near  PI  in  figure  5d  are  also  associated  with  high  radon  levels;  how¬ 
ever,  while  these  trajectory  paths  may  at  first  appear  similar  to  those  corresponding  to  the  low  radon  category  in 
figure  4,  they  are,  in  fact,  distinctly  different  in  that  they  turn  into  the  southern  California  bight  just  to  the  north  of 
San  Miguel  Island.  It  is  noteworthy  that  the  two  northwesterly  trajectories  observed  in  figure  6d,  during  the  tran¬ 
sition  regime,  correspond  to  radon  levels  of  14.1  and  14.7  pCi/m3  while  the  other  two  trajectories,  which  traverse 
the  bight  region,  correspond  to  radon  values  of  25.3  and  27  pCi/m3.  The  two  westerly  trajectories  shown  in  figure  6e 
ended  at  SNI  simultaneous  with  radon  measurements  of  10.5  and  12.5  pCi/m3. 


Table  2.  Trajectory /Radon  Sample  Sire  for  Radon  Concentration  Categoric* 
During  Each  of  the  Five  Experiment*. 


Measurement  Period 

Radon  Concentration  (pCi/m3) 

<10 

>10,  <30 

>30 

780508-780522 

11 

44 

24 

781106-781117 

1 

6 

9 

790425-790511 

54 

15 

1 

800113-800125 

16 

4 

41 

800721-800731 

45 

2 

0 

TOTAL  SAMPLE  SIZE 

127 

71 

75  =  273 
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Figure  4.  Air  trajectories  terminating  at  SNI  during  each 
measurement  experiment  are  illustrated  for  the  low 
(<10  pCi/m 3)  radon  category.  In  general,  low  radon 
levels  measured  at  SNI  are  associated  with  a  mostly 
overwater  northwesterly  to  westerly  trajectory  pattern. 
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Figure  5(a). 
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Figure  5(b). 
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Figure  Sic). 
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Figure  5(d). 


Figure  S.  Air  trajectories  ending  at  SNI  during  each  measurement  experiment  are  illustrated  for  the  high 
(>30  pCi/m3)  radon  category.  High  radon  levels  were  not  observed  during  the  last  experiment  (July 
1980).  In  general,  high  radon  levels  measured  at  SNI  are  associated  with  air  trajectories  that  have  a 
history  over  land  or  in  the  southern  California  bight. 
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Figure  6(b). 
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Figure  6.  Air  trajectories  terminating  at  SNI  during  each 
measurement  experiment  are  illustrated  for  the  medium 
f>10,  <30 )  radon  category.  These  maps  exhibit 
trajectory  paths  that  are  simitar  to  those  found  in  both 
the  low  and  high  radon  concentration  categories  and  are 
suggestive  of  a  transition  mode. 
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Figure  6(e). 


S.  SYNOPTIC  CONDITIONS  AND  RADON  LEVELS 


The  most  striking  difference  in  trajectory  paths  occurs  between  the  low  (figure  3b)  and  high  radon  (figure  3d) 
concentration  categories.  Some  of  the  essential  and  representative  meteorological  features  responsible  for  this  bifur¬ 
cation  will  now  be  discussed.  First  we  attempt  an  explanation  of  the  two  “anomalous”  cases,  i.e.,  the  events  where 
relatively  high  radon  levels  were  associated  with  overwater  trajectories,  by  recourse  to  the  meteorological  conditions 
that  produced  them. 

a.  May  1978 

In  all  of  the  periods  studied,  good  agreement  was  observed  between  trajectory  paths  and  radon  measurements 
except  during  the  May  1978  experiment  when  high  radon  levels  were  recorded  for  two  episodes  of  overwater  trajec¬ 
tories.  The  first  episode,  which  consists  of  the  southwesterly  trajectory  paths  exhibited  in  figure  5a,  occurred  on 
May  14,  1978  when  a  well  defined  “Catalina  eddy,”  illustrated  in  figure  7,  became  established  in  the  southern  Cali¬ 
fornia  bight  in  the  wake  of  a  “Santa  Ana”  (offshore  flow)  condition.  A  number  of  factors  could  have  contributed  to 
the  measurement  of  high  radon  levels  coincident  with  the  oceanic  trajectory  paths.  Four  possible  explanations 
follow:  (1)  radon  may  have  been  advected  from  land  and/or  from  the  post  “Santa  Ana”  air  by  the  eddy  circulation 
while  limited  wind  information  over  the  ocean  could  have  resulted  in  inaccurate  trajectories; (2)  above  background 
radon  exhalation  from  southwest  SNI  could  have  been  transported  to  the  measurement  site  by  the  observed  trajec¬ 
tories;  (3)  if  the  trajectories  were  extended  backwards  far  enough  they  may  reveal  an  overland  history  that  could 
account  for  all  or  part  of  the  high  radon  values  (the  trajectories  terminate  near  the  southwest  map  boundary  when 
extended  out  to  48  hours);  (4)  the  observed  surface  air  trajectories  may  be  accurate  but  not  completely  representa¬ 
tive  of  the  air  mass  history  due  to  horizontal  and  vertical  mixing  effects  that  have  been  neglected  by  our  treatment. 
In  any  event  it  appears  likely  that  the  cyclonic  (counterclockwise)  eddy  circulation  transported  radon  from  land 
and/or  from  post  “Santa  Ana”  air  and  through  the  bight. 

The  second  episode  during  which  high  radon  activity  was  associated  with  overwater  trajectories  occurred  about 
one  week  later  and  is  distinguished  by  the  two  northwesterly  paths  in  figure  5a.  The  pacific  high  pressure  system 
west  of  Washington,  shown  in  figure  8a,  apparently  circulated  air  along  an  overland  coast  route  where  high  radon 
levels  could  be  accumulated,  then  turned  offshore  tracing  out  the  northwesterly  path  illustrated  in  figure  5a. 

b.  November  1978 

The  trajectory  paths  for  the  November  1978  measurement  period  (figures  4b,  5b,  and  6b)  are  somewhat  simi¬ 
lar  to  each  other  and  were  associated  with  high  pressure  off  the  coast  extending  inland  over  northern  California  and 
Oregon  (figure  8b).  Offshore  pressure  gradients,  and  therefore  higher  radon  concentrations,  were  established  due  to 
inland  ridging  to  the  north  in  the  wake  of  transient  short  wave  troughs.  The  air  traversing  the  northern  coast  still 
retained  some  continental  (higher  radon)  characteristics  after  its  short  over-water  transit  to  SNI. 

c.  April/May  1979 

Figure  5c  shows  only  a  single  trajectory  for  the  high  radon  concentration  category  during  the  April/May  1979 
experiment.  The  synoptic  features  producing  this  pattern  are  illustrated  in  figure  8c  and  are  similar  to  those  dis¬ 
cussed  in  section  5b.  A  weak  offshore  pressure  gradient  developed  over  southern  California  as  the  semipermanent 
Pacific  high  was  distorted  by  the  short  wave  passage  over  the  ridge. 

Most  of  this  period  was  dominated  by  the  strong  northwesterly  trajectory  pattern  shown  in  figure  4c.  These 
trajectories  were  supported  by  an  intense  high  pressure  system  located  over  the  Pacific.  Figure  8d  shows  an  example 
of  the  synoptic  conditions  when  radon  levels  at  SNI  were  2  pCi/m3  and  the  trajectories  from  the  regime  in  figure  4c. 
The  lowest  radon  levels  for  all  five  measurement  periods  were  observed  during  this  episode. 


Figure  80.  24  Jui  80  11200  2) 
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Figures  4d  and  5d  illustrate  the  greatest  dichotomy  in  trajectory  paths  for  an  individual  measurement  period. 
The  low  radon  trajectories  in  figure  4d  were  associated  with  the  movement  of  a  frontal  system  across  the  northern 
fringe  of  the  Pacific  High  as  shown  in  figure  8e.  As  the  front  moved  southeastward  and  passed  through  the  southern 
California  area,  onshore  pressure  gradients  with  northwesterly  trajectories  and  low  radon  levels  were  maintained. 

The  high  radon  trajectories  in  figure  5d  resulted,  for  the  most  part,  from  a  strong,  high  pressure,  build  up  over 
the  Great  Basin  that  established  weak  offshore  pressure  gradients  over  southern  California.  This  weather  system  is 
outlined  in  figure  8f  and  was  responsible  for  relatively  stagnant  air  circulation  in  the  bight. 

e.  July  1980 

Only  low  radon  measurements  with  westerly  trajectory  paths,  shown  in  figure  4e,  were  observed  during  the 
July  1980  experiment.  The  air  mass  conditions  for  this  episode  were  characterized  by  air  arriving  from  the  periphery 
of  the  Pacific  High  blowing  towards  the  thermal  trough  over  inland  southern  California  as  illustrated  in  figure  8g. 
Weak  onshore  pressure  gradients,  sustained  by  the  combination  of  the  Pacific  High  and  thermal  trough,  produced  the 
light  westerly  winds  and  low  radon  levels  measured  at  SN1. 


6.  RADON  VERSUS  INVERSION  ALTITUDE 

Rawinsonde  data  from  SNI  was  used  to  construct  time  series  plots  of  potential  temperature  for  altitudes  from 
the  surface  up  to  6000  ft.  The  base  of  the  marine  inversion  layer  was  then  approximated  as  the  height  at  which  isen- 
tropic  packing  was  observed  on  the  time  series  graphs.  Figure  9a  shows  the  scatter  diagram  obtained  for  radon  con¬ 
centration  as  a  function  of  inversion  altitude  for  all  measurement  periods  when  data  was  available.  Evidently,  a 
general  trend  for  higher  radon  levels  is  indicated  when  the  height  of  the  inversion  base  is  low,  and  only  low  radon 
values  are  observed  for  the  largest  inversion  altitudes.  Furthermore,  figure  9a  also  indicates  a  lack  of  dependence 
between  radon  activity  and  inversion  height  when  radon  levels  are  less  than  10  pCi/m3.  The  data  stratum  in  low 
radon  activity  and  the  general  clustering  of  data  in  figure  9a  are  suggestive  of  an  interdependence  between  radon  and 
inversion  altitude  and  warrants  further  attention. 

Figures  9b  and  9c  illustrate  the  trajectory  paths  when  the  inversion  altitude  is  above  750  ft.  and  radon  levels 
are  <10  pCi/m3  and  >10  pCi/m3  respectively.  Figures  9d  and  9e  show  the  trajectories  for  the  same  radon  categories 
but  for  inversion  altitudes  below  750  ft.  From  figures  9b  and  9c,  we  generally  find  high  inversion  (>750  ft)  condi¬ 
tions  corresponding  to  swift  northwesterly  trajectories  with  low  (figure  9b)  and  high  (figure  9c)  radon  levels  depend¬ 
ing  upon  the  trajectory  route.  In  contrast,  low  inversion  conditions  are  mostly  associated  with  high  radon  concentra¬ 
tions  and  upstream  stagnation  in  the  bight  (figure  9e);  a  few  cases  of  low  radon  overwater  paths  also  occur  (figure 
9d).  Since  the  distribution  of  atmospheric  radon  generally  decreases  exponentially  with  altitude  (Jacobi  and  Andre. 
196346;  Larson  and  Hoppel,  197347),  the  higher  inversion  altitudes  may,  in  some  cases,  allow  radon-rich  air  near  the 
land  surface  to  mix  with  radon-poor  air  at  higher  elevations.  Low  inversion  conditions,  by  comparison,  favor  the  sup¬ 
pression  of  deep  vertical  mixing  and  provide,  in  the  presence  of  weak  offshore  pressure  gradients,  ideal  conditions 
for  high  radon  levels  and  the  relatively  stagnant  trajectories  observed  in  the  southern  California  bight  (figures  5a.  5c, 
5d).  In  order  to  completely  account  for  mixing  processes,  3-dimensional  air  motions  over  longer  time  scales  should 
be  considered. 


4^ Jacobi.  W.,  and  K.  Andre,  1963:  The  vertical  distribution  of  Radon  222,  Radon  220  and  their  decay  products  in  the  atmosphere. 
J.  Geophys.  Res.,  68,  3799-3814.  UNC  LASSIFIED 

47Larson,  R.  E„  and  W.  A.  Hoppel,  1973:  Radon  222  measurements  below  4  Km  as  related  to  atmospheric  convection.  Pageoph., 
105.  900-906.  UNC  LASSII  II  D 
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Figure  9(d).  Figure  9(e). 

Figure  ft  Scatter  diagram  thowt  relationship  between  radon  concentration  and  the  approximate  altitude  of  the  bate  of  the 
marine  inversion  layer  for  available  data  during  the  five  experiments.  Figures  9<bl-9(al  illustrate  trajectory  paths  isolated 

according  to  radon  category  and  inversion  altitude. 
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7.  RADON  VERSUS  AIR  PARCEL  SPEED 


Some  of  the  24-hour  trajectories  that  are  coincident  with  low  radon  activity  terminate  off  the  map  boundaries 
(figure  9b),  indicating  high  wind  speeds,  while  other  paths  that  are  associated  with  high  radon  concentrations,  reveal 
much  shorter  histories  (figure  9e)  and  therefore,  lower  wind  speeds.  To  study  this  relationship  further,  we  com¬ 
puted  the  average  air  parcel  speed  along  trajectories  for  up  to  48  hours.  Plots  of  radon  concentration  versus  average 
wind  speed  for  trajectories  of  various  durations  (3,  6,  9  ...  48  hrs.)  show  a  tendency  for  high  radon  levels  to  occur 
at  low  wind  speeds  and  low  radon  levels  to  occur  at  high  wind  speeds.  Figure  10a  shows  the  results  for  1 2-hour  tra¬ 
jectories  during  all  measurement  episodes  and  figure  10b  illustrates  the  contribution  of  the  April/May  1979  period 
to  figure  10a.  Each  measurement  period  exhibited  the  same  trend  except  for  the  July  1980  experiment  (figures  4e, 
6e)  when  low  radon  levels  were  observed  with  light  westerly  winds. 


8.  PSEUDO  MEAN  RADON  DISTRIBUTION 

Figure  1 1  was  produced  in  an  effort  to  parameterize  some  of  the  relationships  that  we  have  observed  between 
air  trajectories  and  radon  measurements  at  SN1.  This  map  illustrates  the  distribution  of  average  radon  levels  for  1 2- 
hour  trajectories  passing  over  the  regions  where  the  values  are  indicated.  For  example,  figure  1 1  shows  us  that  the 
average  radon  concentration  measured  at  SN1  for  trajectories  that  pass  over  Santa  Catalina  Island  (see  figure  1) 
before  arriving  at  SNI  is  approximately  50  pCi/m3.  As  observed  earlier,  but  now  more  quantitatively,  we  find  the 
highest  average  radon  levels  measured  at  SNI  occur  for  air  parcels  passing  over  land  or  through  the  bight;  whereas  the 
lowest  levels  occur  for  air  that  passes  over  water  from  the  west  and  northwest  of  SNI. 

Because  any  air  parcel  traced  sufficiently  backward  in  time  will  ultimately  traverse  the  open  ocean,  the  tech¬ 
nique  discussed  here  should  be  limited  to  short  air  parcel  histories.  The  choice  of  the  12-hour  trajectories  was  a 
compromise  based  upon  this  constraint.  Trajectories  of  24-  and  48-hour  durations  reveal  the  same  trend  in  average 
radon  distribution,  but  with  different  magnitudes.  In  the  future  we  will  attempt  to  quantify  these  results  further  by 
deriving  effective  radon  emission  rates  in  lieu  of  the  average  radon  values  in  figure  1 1 ,  based  on  the  distribution  of 
trajectory  paths  and  the  single  station  radon  measurements. 


9.  CONCLUSIONS 

Despite  the  simplicity  of  the  trajectory  technique  used  here,  which  neglects  topographical,  vertical/horizontal 
mixing,  and  dynamical  effects,  and  considering  the  uncertainties  involved  in  reconstructing  air  parcel  history  from  the 
scarce  wind  information  available,  the  trajectories  are  in  general  agreement  with  the  radon  measurements.  The  bulk 
of  the  radon  fluctuations  are  explained  by  trajectory,  although  inverse  correlations  between  radon  and  inversion 
height,  and  radon  and  wind  speed  were  also  observed.  Because  sufficient  wind  data  for  constructing  air  trajectories  is 
not  available  in  many  areas,  radon  is  a  valuable  diagnostic  tool  that  can  be  used  in  realtime  operational  situations  or 
in  subsequent  analyses.  However,  if  a  predictive  capability  is  desired  for  assessment  of  electro-optical  conditions,  air 
mass  trajectories  constructed  from  forecast  wind  fields  must  be  employed.  We  believe  that  some  effort  should  be 
devoted  to  developing  an  air  mass  index,  based  upon  trajectory,  which  can  be  used  in  aerosol  models  for  predicting 
electro-optical  extinction  in  the  same  fashion  that  radon  is  used  as  a  realtime  air  mass  parameter  (Gathman,  198219). 
Our  results  suggest  the  use  of  marine  aerosol  models  for  predicting  extinction  coefficients  when  radon  levels  are  less 
than  10  pCi/m3;  for  radon  levels  ^30  pCi/m3,  continental  aerosol  models  would  be  appropriate.  Caution  should  be 
exercised  here  since  these  results  apply  only  to  SNI  and  should  not  be  used  for  other  locations,  where  geological  fea¬ 
tures  and  crustal  radium  deposits  may  differ  substantially  (Larson  et  al.,  198248),  unless  confirmed  by  some  trajec¬ 
tory  technique.  Similar  trajectory/radon  relationships  can  be  determined  for  other  areas  of  interest  and.  in  this 
fashion,  radon  levels  can  be  “calibrated”  to  the  local  marine/continental  conditions.  To  test  our  findings  and  to 
more  fully  determine  the  impact  air  mass  effects  have  on  aerosol  properties,  future  investigations  should  explore  the 
feasibility  of  using  radon  and  air  trajectories  to  characterize  aerosol  size  distributions. 


4 8 Larson,  R.  E„  D.  J.  Bressan,  and  P.  E.  Wilkniss,  1982:  Atmospheric  measurements  by  the  Atmospheric  Physics  and  Chemical 
Oceanography  Branches.  NRL  report  in  press,  1982.  UNCLASSIHF.D 


Figure  1 1.  Plotted  numbers  and  isopleths  indicate  average  radon 
concentrations  (measured  at  SNI)  for  12-hour  trajectories  (ending  at 
SNI)  passing  over  the  indicated  areas. 
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